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Abstract. Transient Elastography enables detection and characterization 
of tissue abnormalities. In this paper we assume that the displacements are 
modeled by linear isotropic elasticity system and the tissue displacement has 
been obtained by the first step in hybrid methods. Then we reconstruct the 
Lame parameters of this system from knowledge of tissue displacement. We 
show that for a sufficiently large number of solutions of the elasticity system 
and for an open set of the well-chosen boundary conditions, (A, fi) can be 
uniquely and stably reconstructed. The set of well-chosen boundary con- 
ditions is characterized in terms of appropriate complex geometrical optics 
solutions. 

1. Introduction 

Medical imaging is the technique and process used to create images of the 
human body for clinical purposes or medical science. Each available imagine 
method has its advantages and disadvantages. For instance, in breast imaging 
ultrasound provides high resolution, while suffers from a low contrast. Medi- 
cal imaging modalities such as Computerized Tomography (CT) and Magnetic 
Resonance Imaging (MRI) are examples of modalities providing high resolu- 
tion. In some situations, these modalities fail to exhibit a sufficient contrast 
between different types of tissues. Other modalities, based on optical, elastic, 
or electrical properties of these tissues, display high contrast, such as Optical 
Tomography (OT) and Electrical Impedance Tomography (EIT). For example, 
many tumors absorb more energy of electromagnetic waves than healthy cells. 
This means that using such electromagnetic waves offers high contrast, but 
the resolution is low in this case. 

In order to obtain better image, the natural idea is to try to combine dif- 
ferent imaging modalities. These are coupled-physics imaging methods, also 
called hybrid methods. It is to combine the high resolution modality with 
another high contrast modality. Examples of possible physical couplings in- 
clude: Photo-Acoustic Tomography (PAT) and Thermo- Acoustic Tomography 
(TAT), Ultrasound modulated Optical Tomography (UMOT) and elasticity 
with ultrasound in Transient Elastography (TE). 

Reconstructions in hybrid inverse problems involve two steps. The first 
step is to solve the inverse problem concerning the high-resolution-low-contrast 
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modality. For instance, in PAT and TAT, this corresponds to the reconstruc- 
tion of the map of absorbed radiation from the ultrasonic measurements from 
boundary acoustic wave measurements. In Transient Elastography, this is to 
solving an inverse scattering problem in a time-dependent wave equation. In 
the second step, we assume that the first step has been done and we try to 
reconstruct the coefficients that display high contrast from the mappings ob- 
tained during the first step. 

In this paper, the modality we consider is Transient Elastography. See [10] 
and references there for more details. In TE, the high resolution modality is 
also ultrasound. The tissue displacement data can be obtained by ultrasound 
in the first step. The second step is to recover the Lame parameters from the 
knowledge of the tissue displacement. In the following paper, we will assume 
the first step has been performed. 

The displacement is modeled by the following equations. Let Q C M n , n = 
2, 3, be an open bounded domain with smooth boundary. Let u be the dis- 
placement satisfying the linear isotropic elasticity system 



where S(A) = (A + A T )/2 is the symmetric part of the matrix A. Here (A, jj) 
are Lame parameters and w 6 Mis the frequency. The set of internal functions 
obtained by ultrasound in TE is given by : 



Denote that V = {(A,/i) G C 7 (Q) x C 9 (Q); < m < ||A|| c7( n), IMI^p < 
M and A, > 0}. Let H and H be two sets of internal data of the elasticity 
system with parameters (A,/i) and (A,/2), respectively. Below is our main 
result. 

Theorem 1.1. Let Q be an open bounded domain ofMJ 1 with smooth boundary. 
Suppose that the Lame parameters (A,/i) and (A,/t) G V and fi\gn = fi\dn- 
Then there is an open set of the given boundary data such that the uniqueness 
and the global stability of (A, /i) hold. 

The remainder of this paper is organized as follows. In section 2, we in- 
troduce the CGO solutions of the elasticity system. Section 3 is devoted to 
constructing the Lame parameters in two-dimensional case. The reconstruc- 
tion of (A,/i) in 3D is presented in section 4. 

2. Complex geometric optics solutions of the elasticity system 

In this section, we will briefly introduce the complex geometric optics (CGO) 
solutions of the elasticity system. Based on [9], we can derive the following 
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Here Vq contains the third derivative of \i and 

-2/i 1 / 2 W 1 + /i" 3 ^ 2 - V log fi 

Then the solution of the elasticity system (1) is 



VAx) 



u ■= [T 1/2 w + /i _1 V/ - fVfi~\ 

Here V 2 g denotes the Hessian matrix d 2 g / dxjdxk- Note that we will not need 
the explicit form of Vq in the construction of CGO solution. 

The construction of CGO solutions of (2) with linear phase was first deduced 
by Nakamura and Uhlmann in [11] and [12], where they introduced the inter- 
twining property in handling the first order term. Eskin and Ralston [7] also 
gave a similar result in 2002. Later Uhlmann and Wang [15] used Carleman 
estimate to deduce the CGO solutions in two-dimensional case. 

In the following sections, we will use the CGO solutions constructed in Eskin 
and Ralston's paper [7]. For the rest of this section, we will introduce the key 
lemma from [7]. 

To construct CGO solutions of (2), it is convenient to work on M n instead 
of Q. Since Q is bounded, we pick a ball Br for R > 0, such that Q C Br and 
extend A and fi to R n by preserving its smoothness and also suppA, supp/x C 

Br. 

Let a and /3 be two orthogonal unit vectors in M. n . Denote that p = r(a + 
ij3) G C" and 9 = a + ij3 with r > 0. The following lemma proved by Eskin and 
Ralston in [7] in three-dimensional case, but it is still true in two dimensional 
space. 

Lemma 2.1. (Eskin- Ralston) Consider the Schrddinger equation with external 
Yang-Mills potentials 

(3) Lu = -Au - 2iA(x) ■ Vu + B(x)u = 0, x G B R C lR n 

where A(x) = (Ai(x), . . . ,A n (x)) G C no (B R ), n > 6 with Ai(x) and B(x) are 
{n + 1) x (n + 1) matrices. Then there are solutions of (3) of the form 

u = e ip - x (C (x, 9)g{9 ■ x) + 0(t~ 1 )) 

where Cq G C n °(BR) is the solution of 

i9-^-C (x,9)=9-A(x)C (x,9) 

with det Cq is never zero, g(z) is an arbitrary polynomial in complex variables 
z, and 0(t~ 1 ) is bounded by C(l/r) in H k (B R ), < k < n - 2. 

By the lemma above, the CGO solutions of (2) can be written as follows 

W ^ ( r ) +0(r~ v 



f J \\s 

where (r, s) T is C (x,9)g(9 ■ x). We can write w = e lp ' x {r + 0(r -1 )) and 
/ = e i P- x {s + 0{T- 1 )). Plug (w,f) T into u = yT l / 2 w + f - fV yT l . Then 



4 



we have the CGO solutions of the elasticity system. Note that (r, s) T satisfies 
the equation 

(4) -29 ■ V ([ ) = Vy{x) 
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Remark 1. Since Cq is invertible at every point in Q, we can conclude that 
(r, s) T is not zero everywhere in Q. This does not imply that both r and s 
never vanish in Q. However, for any point y e Q, there is a small neighborhood 
By of y in fl and a CGO solution of (w, f) T such that both r and s do not 
vanish in B y . 

3. Reconstruction of Lame parameters in two-dimensional case 

In the previous section, we already have the CGO solutions of the elasticity 
system. Now we want to use them to give a reconstruction of \i. Let u = 
(ui,U2) T be the displacement which satisfies the elasticity system 

(5) V- (A(V-u)/ + 2S(Vu)/i)+w 2 u = 0. 

From (5), we can deduce the following equation 

fd l (V-u)\ I A + ,u \ ( a + b \ 
A + /i 
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di(X + fi) di(a 
\ 9 2 (A + /i) J \ d 2 (a 

where a = d^Ux + dxU2, b = dxUx — d 2 U2, and u* = {ux + u 2 ). 

Inspired by the paper of Bal and Uhlmann in 2012 [3], we can eliminate the 
first term of (6). Since the vector 

/ dx(V-u) \ 
d 2 (V-u) 

V • u 

V • u 



(7) 



V 



just has three different entries, we only need to construct three linearly inde- 
pendent vectors on some subdomain in Q. With these three vectors, we can 
remove the first term of the left hand side of the equation (6). More precisely, 
suppose that vP, for j = 1, 2, 3, are three different solutions of (5). Then 



/ 0i (V 
d 2 (V 



satisfies that 
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A + /i 
A + fi 

\ d 2 (X + fi) J 
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for j = 1, 2, 3, where u J * = u{ + u 2 . Assume that {v\, v 2 , ^3) are three linearly 
independent vectors in some subdomain of Q, say Q,q. Then there exist three 
functions 61, 62, and 63 such that v can be expressed as v + Y^=i ®j v j = 0- 
Multiplying (8) by Qj and summing over j with equation (6), we have 



3°i> 



d 2 p 
fi 
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Let /3 = (v • ei)ei + (v • e 2 )e 2 and 7 = (e 3 + e 4 ) • v. Here v = v + X/j = i ©j^j and 



e,- G 



with the j entry is 1 and others are zero. Then the above 



equation can be rewritten as 



(5 ■ V/i + 7/i 



2/ > 
-CO [U 



3=1 



Suppose that the real and imaginary part of /3(x) are linearly independent in 
R 2 for every x G Qo . Then we can recover p in Q C for each frequency u; 
independently. 

3.1. Local stability of p in 2D. We will show that how we can get three 
linear independent vectors of the form (7). We plug the CGO solutions u p = 



P 



-1/2 



w P + n l Vf p - fpVff 



1 into v p , where 
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Then we get v p is 
(10) 
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-1/2/ 



/ <9i(V-u p ) \ 

V • u p 

V • U n 



2fi l {p\d 1 s + pip 2 d 2 s) 



-p- 1/2 (p 2 2 r 2 + pifrn) - 2p~ 1 {p 2 2 d 2 s + Pl p 2 d lS ) + 0(\p\) 



p 2 r 2 ) + 2ip 1 p ■ Vs 
ip^ 1/2 (p 1 r 1 + p 2 r 2 ) + 2z'/i _1 p • Vs 



+ 0(1 
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Since (A, /i) G "P, by Lemma 2.1, it implies that r and s are in C 7 and C^r" 1 ) G 
H 5 . Note that r = (r 1; r 2 ) T and p = (pi,p 2 ) T . 

Now we fix any point x G <9fi. Let p = r(l, i) = — ip G C 2 with r > 0. 

Since, in Lemma 2.1, the matrix solution Cq(x, 9) is invertible, we can choose 
a constant vector such that Cq(x, 9)g(°' = (r^°\ s(°)) T with the conditions 

s (0) (x ) = 1, s {0) ^ and p • r {0) ^ 

in a neighborhood of £o m Br, say {7. Then we have the CGO solution of the 
elasticity system, that is, 



u 



(0) 



/^v/f - /j 0) Vp 



6 

where 
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and C2{x,6) satisfy the following two 

i9 ■ ^C 2 (x, 8) = 8 ■ V^C^x, 8), 

respectively. Since p = ip, we can choose C2{x,8) = Ci(x,9). Moreover, with 
suitable constant vector g, we can get that r^> = and = and 

sW(x ) = , \VsW(x Q )\ < y/2/8; r«(x ) = (1, -i) = r^(x ). 

By continuity of r^ k \ we have p ■ ^ in a neighborhood U of Xq. Then the 
CGO solutions are 



Let 8 = p/r,6 = p/r. Let d(x,6 
equations 

id~C 1 (x,0) = 6-V 1 (x)C 1 (x,6), 



u 



(fc) 



where 
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r (l) 

s (l) 



+ 0(r 



(2) 
P 

(2) 



r (2) 
s (2) 



)+0(r- 1 )). 



So far we have three CGO solutions, that is, u p °\ u p l \ and u p 



(2) 



From the 



previous section, the equation (4) gives us 
(11) 



2p . Vs = ^V /2 p • r. 



X + 2p 

Using (10) and (11), we obtain that 



u p := e 



-ip-x 



\ P \-y 



\+2p 1 
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(i) Jih 



0(i 



where = (r^ , ). Moreover, at point xo, 

/ 2 \ 



uJ(a:o) 



1 v 7 / 1 



2A + 2/x 



2i 


V / 
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Similarly, we use (10) and (11), then we obtain that 



/ _^L r ( r (D +i r W) + 0(1) \ 



u p := e 



(i) Jih 



X+2p 



0(1) 



l X+2 IJ ,\ r i 



ir 



(i)> 



0{t 



and 



V2A + 2/i 



( -2r + 0(l) \ 
-2ri + 0{l) 
2i 
2i 



Oir- 1 ) 



Also, 



u 



and 



V2 



/ ^/ _(2) , (2). 

_ VP T (-j r {2) +r {2h 



Z A+2^I^ r l r 2 J 

; /, v (2) _ _(2)n 



0(1) \ 
0(1) 

) 



Up(^o) 



V/ 1 



V2A + 2/1 



/ 2r+0(l) \ 
2ri + 0(l) 
-2 
-2 



oCr- 1 ; 



Let u p = u 2 , + Up. Since r^ 2 ^ = and s^> = we have 



u, 



V2 



; VP 
' A+2p 



0(1) 
0(1) 



((l + z)ri 1} + (-l + z)r^) 



+ 0(0 



and 



1 Vg 
A + 2/x 



/ 0(1) \ 
0(1) 
2i - 2 

y 2? - 2 y 



+ 0(0. 



We denote 3ft/ and S/ to be the real and imaginary part of /, respectively. 
Then 9ftu p , ^sut, and 9ftu p are linearly independent in the neighborhood U 
of Xo when r is sufficiently large. Therefore, there exist functions 9j and 
j = 1, 2, 3, such that 



9fte-^°/|p| 2 + ^!3ftu p + # 2 3u p + ^ 3 3ftu„ = 



and 



$Se- ip - x v° p /\p\* + «u p + tf 2 5u* + tf 3 9forf = 0, 



respectively. 

Summing over j, we have the following equations 



Pi ■ V/i + 7i/i = -uj 2 Ui 



for z = 1,2, 
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where 



IpI 2 V a - b J IpI 2 V a b J IpI 2 \ a b 

e~ ip - x ( a 1 + 6 1 \ e-^ x f a 2 + b 2 



\p\ V al - &1 ; ~W \^ 2 - b2 



and 

7i = 



e-^' x / a 1 + 6 1 \ e"^ x f a 2 + b 2 



IpI V a - 6 / IpI V a 6 

72 = ( j ) • ( g ) + + ^ + • 



Note that 



i _ _ ip . X| ,_ 2 / ^(a 1 + 6 r 



and 



Moreover, 



-ip-x p -ip-x p -tp-x 

IpI 2 IpI IpI 

p—ip-x p—ip-x 

eM—-u^ + e — r u^ 
\p\ \p\ 



and 



p — lp-X p — ip-x p — ip-x 

IpI 2 IpI 2 IpI 2 

(p-ip-x p-ip-x 
IpI IpI 

where u^ 7 ')* = + ?4 . 

Remark 2. By choosing suitable |03(#o)| an d I^O^o)! can be less than 
or equal to 1. To show that, we let m = max{|6 l 3(a;o)|, |$3(xo)|}- We choose a 
new constant vector g(°\ instead of the orig inal £> (0) , such that (f^, s (0) ) T = 
C (x,6)gW(6-x) where r^(x ) = r<®{x )/M, M > m, §(°\xo) = s^(xq) and 



(f'°',^°') T satisfies the original assumption, that is, s^°\xo) = 1 and ^ 
and p ■ f ^ in a neighborhood of Xq. Then 

■3^(1.0 -^ 0) (xo) 



to) 



-P 3j 3ft 



/ / 

1 
2 



V 



V 



o 
o 



+ o(r- 1 : 



/ 



Pi 



/ 



and 



§ 3 (x ) = -p 3 m 



( ( 

i 



2M 



V 



■^(M)-r<°>(*o) \ 



V 








oo-- 1 ; 



Pv 



J 



where P = (pif) with pij = p^ ■ pj is a boundedly invertible symmetric matrix 
with P u the entries of P" 1 . Here p\ = Sftu*, = Su*, ^3 = 9ftu^. From 

above, we obtain that #3 is close to 63/ M as r is large. Therefore, the new 
|#3(xo)| are less than or equal to 1 when r is sufficiently large. Similarly, we 
have |t?3(xo)| < 1. In fact, |^3(xo)| and |$3(#o)| can be very small when we 
choose a suitable constant vector g(°\ 

Lemma 3.1. (3j is not zero in some neighborhood of Xq, for j = 1,2. 

Note that we denote by ~ equalities up to terms which are asymptotically 
negligible as r goes to infinity. 

Proof. From above construction, we have knew that 6^ ~ r° and i?j ~ r°. At 
point x , since s^(x ) = 0, we have 



(12) 
and 
(13) 



9ft 



3 



-ip-x 



-tp-x 



\p\ 



b 1 
b 1 



b 1 
b 1 



[X ~ T 



Since s^ k \x ) = 0, r^ = = and (1, —i) ■ r^(x ) = 0, we obtain 

that 



9ft 



3ft 



-ip-x 



a 1 + b 1 
a l -b l 



-lp-X 



a 2 + b 2 
a 2 -b 2 



2^ 1 r(2id 1 s^ + 2<9 2 s (1) ) 
2^ 1 r(-2d 1 s^ + 2id 2 sM) 



V2t 

4 „-i» 1 /7 + 



(a: ) 



+ 0(1) (so) 



0(i 



^0 
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Under the conditions |Vs (1) (x )| < y/2/8, \9 3 (x )\ < 1, and s (0) (x ) = 1, then 
we have 




v* / a 1 + b 1 



\f> 



^(_d lS (D + ^ 2S (i) 




Combining (12), (13) and the asymptotic result above, we have 




lifts* 1 * + 




By continuity of the result follows. □ 

Based on the above lemma, we may suppose that (3j is not zero in the 
neighborhood of xq, say V. Moreover, with the suitable chosen CGO solution 
u, these two functions 9j and could be small enough such that f3\ and /3 2 are 
linearly independent in 7 as r sufficiently large. Then we have the following 
equation 

(14) Pj ■ V/i + ij/i = j = 1, 2, 

with {f3j(x)}j=i t 2 a basis in M 2 for every point x G V. 

Thus, there exists an invertible matrix A = (a^) such that = a^efe 
with inverse of class C°(fl). Thus, we have constructed two vector- valued 
functions T(x), $(x) G (C°(f2)) 2 such that (14) can be rewritten as 

(15) VyU + rXx)/! = <$>(x) mV. 

Let flo = [/ H 1/ fl O and fT p = (u p , u p , uf 1 ). Then we have the following 
theorem. 

Denote that H p and H p are two sets of internal data of the elasticity system 
with parameters (A,/x) and (A, jl), respectively. 

Theorem 3.2. Suppose that (A,//) and (A,/2) G P. Then, for any fixed xq G 
<9f2, i/iere exists an open neighborhood Qq of Xq in Br such that we have the 
uniquely local reconstruction of \l and the local stability of \i in Q . More 
precisely, we have the estimate 

(16) - Allco(Oo) < C (\fi(x ) - ji(x )\ + \\H P - flp||ca(Oo)) • 



11 

Proof. In the construction of CGO solutions, we take Q C Br. By the argu- 
ment above, the domain Qq can be taken as a small open ball with center xq 
such that Qq C Br. We want to solve the over-determined problem 

j V/i + T(x)fx = in Q , 

\ fi(x ) = /i(x ), x G dtt. 

Let x G flo be an arbitrary point and ip(t) = (1 — t)x + tx, £ G [0, 1], be the 
curve connecting point x and x. Restricted to this curve, then we have the 
uniqueness solution 

^(t)) = ^(V>(0))e-£ r W s)), ^ w * + / ■ * e [0, 1]. 

•/ o 

Similarly, jl satisfies 

f V/i + f(x)/t = <f>(x) in n„, 
\ /i(x ) = /i(x ), x G <9fL 

Restricted to the same curve, it can be expressed as 

/2(V>(<)) =/x(V>(0))e-^ f W s)) ^ w *+ [ ®(ip(s)) ■t/j'(s)ds, t G [0,1]. 

Thus, for all i6fl 0l the difference between /i and /i is 
-/2(z)| <|/x(x ) -M(x )||e-/o r (^)H'(^| 

+ | / $(V>(s)) ■ ^'(s)ds - / $(V>(s)) 

JO JO 

Since fi, jl and T, $ G (C°(f2)) 2 , there is a constant C such that 
Ha* - Allc°(n ) < C (|a*(^o) - At(^o)| + 1 1 -Hp - #Jc 2 (n )) 

□ 

Note that, for any xq G Q, we still have the local stability of /i. The reason 
is that we can still construct three linearly independent vectors near the point 
Xq such that the first term of (6) is eliminated. Then we follow the similar 
idea as in Section 3.1. There exists an open neighborhood B XQ of xq and a 
boundary point Xq G dB XQ . We will get that 

(17) ||/i - Mco { b xo ) < C (|/i(4) - /2(4)l + \\ H P - H p \\cHb xo) ) , 
where B Xo is some open neighborhood of xq. 
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3.2. Global reconstruction of fi in 2D. In the following, we will prove 
that given boundary data, there is an open set of these given boundary data 
such that the stability of \x still holds. That is, the small perturbation of the 
boundary data will not change the stability of [i. 

Let g G H 2 (Q). We consider the following boundary value problem 

V • (A(V • u)I + 2S{Vu)n) + co 2 u = in ft, 
u = g on <9ft. 

Then it is suffices to study the problem 

, , j V ■ (X(V -u)I + 2S(Vu)fi) + uj 2 u = -F in ft, 

^ > \ u = on 9ft, 

where F G L 2 (ft). Let 

(20) B(u,v) = / A(V • u)(V ■ u) + 2[j,S(Vu) : S(Vv) - u 2 u ■ vdx 

Jn 

for w, v G i?o(ft)- Then i? is a continuous sesquilinear form on Banach space 
^(ft). Here A : B = £ 2 =1 £ 2 =1 For u G ^(ft), from (19), we have 

(21) B(u,v)= I A(V ■ u)(V ■ u) + 2(iS{Vu) : 5(Vu) - w 2 u ■ Udx 



(22) 



B(u,v)= [ X(V -u)(V -v) + 2nS{Vu): S(Vv 
Jn 

[ F-vdx, for all v G ^(ft). 



Since A > and /i > and by Korn's inequality, there exists a constant c > 
such that 



/ A(V • «)(V ■ u) + 2jiS(yu) : 5(Vu) - w 2 u • udx 



> / 2fiS(Vu) : S'(Vu) - w 2 n ■ ttda; 
Jn 

>(2cfi - uj 2 )\\u\\ 2 H i {n) . 

If a; is small, then the constant 2c\i — u 2 can be positive. This implies that 
B(u,v) is coercive on if 1 (ft). By Lax-Milgram theorem, there exists unique 
solution u G -f^o(ft) such that 

B(u,v) = (F,v) for all v G if 1 (ft). 

Hence, we obtain the weak solution of (19). Let m be an nonnegative integer. 
Applying regularity theorem, if F G H m (Q), then the weak solution u G Hq(Q) 
is in the space if m+2 (ft). Moreover, we have the estimate 

|p||ijm+2(n) < C||-F||ijro(f2). 

If we let the boundary data g G if 4 (ft), based on the argument above, there 
exists a solution u G if 4 (ft) such that 

V- (A(V • w)/ + 2S(Vw)/i) + w 2 w = in ft, 
u = g on <9ft. 
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By Sobolev embedding theorem, H k (Q) C C m (f2), k > m+n/2, we can deduce 
that both g and u are in C 2 (Q). 

Now, we have CGO solutions Up , j = 0, 1, 2. Let gp = Up \qq be the given 
boundary data. By Lemma 2.1, since (A,yu) G V, we knew that Up G if 4 . Let 
^ if 4 (f2) be the boundary data close to gp , that is, 

yO) _ g^\\ H 4( m) < e for some e > 0, 

then we can find solutions vP> of the elasticity system with boundary data 
gti>. By elliptic regularity theorem and Sobolev embedding theorem, we have 

(23) h^-ufWc^a) <Ce 

for some constant C which is independent of A,//. Then we obtain that 

-v®\\c°{n) < Ce, 
where and Vp are the form of (9). Moreover, for j = 1,4, we have 

and 

pui - 9fUp||c°(n) < Ce. 

Therefore, 5Jui, Sui, and 5ftu4 are also linearly independent when e is suffi- 
ciently small. 

Using the discussion above and Theorem 3.2, we have the following local 
stability result. 

Lemma 3.3. (Local stability of fi) Let Q be an open bounded domain of M 2 
with smooth boundary. Suppose that (A, fi) and (A, fi) G V. Then, for any 
fixed xq G dVt, there exists an open neighborhood Qo of xo in Br and an open 
set of the given boundary data (gp )o<j<2 such that we have the uniquely local 
reconstruction of fi and the local stability of fi in Q . More precisely, we have 
the estimate 

(24) - /2||co(n nn) < C (\fi{x ) - fl(x )\ + \\H - H\\ C 2 { n )^ ■ 

Note that in Lemma 3.3, H = (u {0 \u {1 \u^) and H = (u^°\u^\u^). 
For any point x G fi, Lemma 3.3 still holds in the following sense. 

(25) 

11^ - Mco(b xq ) < C (|/i04) - fx{x£)\ + ||if X0 - H Xq \\c\b xq )) , 4 G dB X0 , 

where B xo some open neighborhood of xq in fi. The two internal data H xo and 
H xo are dependent on the point x . 

The global stability result is stated as follows. 

Theorem 3.4. (Global reconstruction of fi) Let fi be an open bounded domain 
of R 2 with smooth boundary. Suppose that the Lame parameters (A,/x) and 
(A,/i) G V and fign — fiao.- Then there is an open set of the given boundary 
data such that the uniqueness and global stability of \i hold. 
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Note that for the uniqueness of fi, we suppose that the two set of internal 
data are the same, that is, H = H. Since \x is uniquely reconstructed near a 
fixed point xo G dfl under the condition ii\qq = /i|an, from the stability of 11 
in Q, we can obtain that \i = jx in Vt. 

Proof. In section 2, we constructed CGO solutions in a ball Br which con- 
tains f2. First, we consider any point x in dVL. Then we can find an open 
neighborhood B x C B R of x. By Lemma 3.3, we have the estimate 

(26) - p,\\c (B x nn) < C\\H X - H X \\ C 2^ 

since li\ 9 q = fl\ d Q. 

Second, for any point y G tt, since the local stability of Lemma 3.3 is still true 
when xq is in Q, there exists an open neighborhood B y C Q with B y D dQ = 



II A* - AllcW < c (l/ i (l/ + ) - Kv + ) \ + \\ H y - H y \\ C 2(n)J for some y + G dB y . 

Therefore, the compact set f2 is covered by {J x pqB x . Then there exists finitely 
many B x , say, B Xl , . . . , B XN , such that Q C UfeLi B x k - 

Now for arbitrary point z£fl, since fl C UfeLi ^xk> there is an B x . such that 
z G S^.. Suppose that B x . D 9fi 7^ 0, this means that is a neighborhood 
of boundary point Xj. Then, by (26), we have 



such that 



(27) 




For the point x^, since f2 is covered by finitely many subdomain B Xk , after at 
most N — 1 steps, we have 



N 



(30) 




Combing (29) and (30). Then we have 



N 



(31) 




k=l 



Therefore, we have the global stability 



N 




k=i 



□ 
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V • u 

a + b 
a — b J 



<9 2 (A + /i) 
difj, 
d 2 [i 



+ 



V 



/ X + fi \ 
A + /i 

/' 



-w 2 (mi + u 2 ). 



3.3. Global reconstruction of A in 2D. The elasticity system can also be 
written in this form 
(33) 

fV-u\ fd 1 (\ + n)\ f d 1 (V-u)\ 

d 2 (V-u) 
d 1 (a + b) 
d 2 (a-b) J 

As in the reconstruction of p, we will construct CGO solutions such that the 
first term of the equation (33) can be eliminated. 

Plugging the CGO solution u p = p~ l l 2 w p + p~ l V f p — f p Vp~ l into 

/ V • u \ 
V • u 
a + b 
\a-bj 



we get 



lp-X 



+ 0(1) 



ix~ 1/2 p ■ r + 2fi~ 1 p ■ Vs 
pr x l 2 p ■ r + 2/i~V • Vs 
2p~ l (is(p\ + pip 2 ) + cut ■ Vs) + /i _1/2 wi • r 
\ \ 2^\is( Pl p 2 - p\) + u 2 ■ Vs) + pT x l 2 u 2 ■ r ) 

where u x = (pi + p 2 , p\ - p 2 ) and u 2 = (p 2 - pi,pi + p 2 ). 

Let Xq be a point in dfl. Denote that p = r(ei + ie 2 ) and p = ip. 

We choose a constant vector g(°> such that Cq(x, 9)g^ = (r^°\s^) with 



,(o) 



^0, p-r {0) ^0, p-r [[)) {x ) = 1 



(o), 



in a neighborhood of x$, U. Then we get the CGO solution of the elasticity 
system 



u 



(0) 



with 



(0) 



f (0) 



r (0) 
s (0) 



We choose another constant vector such that C±(x, 9)g^ = (r^\ s^) 
with 

s (1) y£ in U and p ■ r^(x ) = 0. 
Then we get the CGO solution of the elasticity system 



u 



(i) 



/ i - 1/ M 1) + ^ 1 v/I 1) -/| 1) v/.- 1 



with 



4 1} 
f (i) 



.(1) 



0(i 
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For k = 2,3, we choose a constant vector such that Ck(x,8)g^ = 
( r ( fe ) ) with 

p ■ r (fc) ^ 

in U. Here we can choose 

r (2) = r (3) s (2) = s (3) 

by taking g^ = g^ and (^(x, #) = Cs(x,9). Then we get the CGO solution 
of the elasticity system 



a 



with 



(2) 
Wp 

(2) 



r (2) 
s (2) 



id 



(3) 
P 

(3) 



r (3) 
s (3) 



+ 0(r 



Then we have that 



— te 



+>*\p\- % vl = 2i»~h 



( \ 



1 + z 
i - 1 



+ 0(0- 



For fc = 2, 3, we have that 



u 



ze -^IpI -1 ^ 
/ ^(l,i). r « 



p ' 



A+2/i 

v 7 ^ 2/i~ 1 r(l + i)^ (2) + z/ 2) 
\ 2//- 1 r(i - l)is^ + z v(2) y 

where = 2/i~ 1 (l + f, 1 - i) ■ Vs^ + /i~ 1/2 (l + i, 1 



i) ■ and 



u 



3 — -ze-^lpl" 1 ^ 



1 



^-l)- (3) 

^-l)-^ (3) 
2 / u- 1 r(-l - + 

(3) -I- i,y(3) 



\ 



/ 



0{r- 1 ] 



\ 2/i 1 r(l — i)is yj; + iv y 

where = 2/i~ 1 (-l+z, z + l)-Vs (3) +^ 1/2 (-l+z, i + l)-r^\ Let = u 2 +u 3 . 
Then the higher order is eliminated. Thus we have 



uj = u 2 p + u 3 



1 

7! 



(1 + z, i — 1) • r 



. r (3) 



^(l + M -l). r (3) 



\ 



2/i~ 1 (2i, 2) ■ Vs^ + /i" 1 / 2 (2i, 2) • r( 3 ) 
\y 2/i~ 1 (-2, 2i) • Vs^ + /i" 1 / 2 ^, 2i) • r( 3 ) ^ 



+ 0(r- 1 ). 



From above, we have constructed three linear independent vectors 9ftu*, Su*, 



and 9fan in U as r is sufficiently large. Therefore, there are three functions 
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9 k , k = 1,2,3, such that Ke^^/M 2 + 01^11* + ^uj + fl-^ = 0. Summing 
over j, we get the following equation 

/ A + ^ \ 

X + P 
P 



(34) 



u> 2 u, 



\ f ) 



where 
u = ft 
and 



p —tp-x p —%p-X p —lp-X 

U (O)* +01 ^1_^ U (1)* + 2 Q«1__ U 



z|p| 5 



ft- 



-ip-x 



%\p\' 



i\p\< 



( d 1 {V-uf)\ 
d 2 {V-u (0) 



i\p\- 



(1) *+# 3 9? 



-ip-x 



i\p\ 



-u 



(2)* 



-ip-x 



i\p\ 



-u 



(3)= 



P , 

o , 6 o> 



9i(a' 
V <9 2 (a°-&°) / 



/ Si(V 



-ip-x 



%\p\' 



it 



/ 9i(V 



-ip-x 



U 



(2)> 
P > 
,(2)> 



(3)> 



<9 2 (V ■ u p 
d 1 (a 2 + b 2 
\ d 2 (a 2 -b 2 ) ) 



-lp-X 



+ 



<9 2 (V ■ u p 3) ) 
d^ + b 3 ) 
V d 2 (a 3 -b 3 ) ) 



k\ = ap — oj 2 u, 



We obtain that 
(35) 
where 

k = (i, i, o, o) ■ (t>; + e^l + e 2 ^ p + e 3 KoJ), 
a = -(i, i, i, i) ■ (»; + ^sroJ + e 2 ^l + £ 3 KoJ). 

Lemma 3.5. k does not vanish in some neighborhood of Xq. 

Proof. Similar argument as Lemma 3.1. From the construction, we have known 
that 9j ~ r° and p ■ r^'(io) = 0. Then 

(36) (1, 1, 0, 0) ■ fliStoJ ~ T~\ (1, 1, 0, 0) • 2 3t>J ~ r" 1 . 

Observe that 



(1,1, 0,0) -9^ 



=3? 



1\ / e-*' 8 / ^(V-nW) 
i; A«>l V<9 2 (V-^) 
The O(r) term of (1, 1, 0, 0) • 9^ is 
1 



d 2 (V-u^) 



(-T^(pf } + P? ] )P {2) ■ r (2) ~ T^(pf } + P2 3) )P (3) • ^ = 



iy/2r A + 2/i vri ' rz /r A + 2p' 

since r^ 2 ) = r^ 3 ^. Thus the leading order term of (1, 1, 0, 0) ■ 3F5t>^ is 0(1), that 
is, 
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where 

i\/2r V 

+ 2ip ■ VrW) + pT 1 {2ip ■ (V(9is«) + V(Sbs (fc) )) + *(pi + p 2 ) As (fe) ) 
+ • Vs (fe) - d 11 ii-H(pP + )s« 

+ 2^ 2 /i-V W • - dnirHiJP + p { 2 k) )s^ 
+ dm-Wiip^V? + + 2zpf V! fc) ) 

+ 9 2 p- 1 / 2(ip w r w + ip W r W + 2ip f) r (* h 

Since 

(m, o, o KW =*(i)-nS^(*(? : «)) = + ° (T_I) 

and (36), we get that 

(37) k(x ) =(1, 1, 0, 0) • (0° + + 6 2 ^l + 3 SffoJ)(z o ) 

(38) =(1, 1, 0, 0) • D°(x ) + (1, 1, 0, 0) • 9 3 ®X>*(x ) + 0(r- x ) 

(39) ._Z^_ + fl8(a?oMa ,). 

Since p E V, <p(xo) is some fixed number. Moreover, we can take small 
such that |6 l 3 (xo)| to be sufficiently small, then we can obtain that 2 (a+2^0 + 
^3(^o)<p(^o) 7^ 0. By continuity of k, there exists a neighborhood V such that 
k never vanishes in V. □ 

Let f2 = U (~)V DO,. We have 3?u^, Shi* and Ku^ are linearly independent 
in O as r is sufficiently large and also k does not vanish in Q . Then it follows 
that 

(40) x = ^-co 2 - inO . 

K K 

In the following theorem, the internal data H p contains all CGO solutions we 
constructed in Section 3. 

Theorem 3.6. Suppose that the Lame parameters (A,p) and (A,p) G V. 
Then, for any fixed xq E dQ, there exists an open neighborhood Qq of xq in Br 
such that we have the local reconstruction of A and the local stability of A in 
Qq. More precisely, we have the estimate 

(41) ||A - A|| c °(n ) < C (\fi(x ) - fl(x )\ + \\H p - #Jc 2 (Q )) • 
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Proof. By theorem 3.2 and equation (40), we have 

||A — \\\c°(n ) <0 ( — p\\c°{n ) + \\H p — #p||c 2 (n ; 

lc 2 (n ; 



<C (\p(x )- p(x )\ + || H p -H p 



□ 



With (40) and the local stability result, Theorem 3.6, we follow the same 
proof as in Section 3.2. Then we can get the following global reconstruction of 
A. 

Theorem 3.7. Let Q be an open bounded domain ofM? with smooth boundary. 
Suppose that the Lame parameters (A,/i) and (A,/i) G V and p\dn — p\dn- 
Then there is an open set of the given boundary data such that the uniqueness 
and the global stability of A hold. 

4. Reconstruction of Lame parameter in three-dimensional case 

In the three-dimensional case, the reconstruction of A and p basically follows 
the same idea as in two-dimensional case. The difference between these two 
cases is we will need more CGO solutions in three dimension in order to get 
four linearly independent vectors locally. 

4.1. Global reconstruction of p in 3D. Let u = (ui,u 2 ,u 3 ) T be the dis- 
placement which satisfies the elasticity system 

(42) V- (X(V -u)I + 2S(Vu)p) + cu 2 u = 0. 

From (42), we can deduce the following equation. 



(43) 



/ &(V-u) 
<9 2 (V • u) 
d 3 {V-u) 

V • u 

V • u 

V ■ u 



\ 



J 



where 6 i7 - = d k u k - diUi 



A + fi 

A + pL 
A + pL 

di(X + fi) 
9 2 (A + /i) 
V d 3 (X + fi) ) 

djUj + diUk 



+ 



^23 
6l3 

bn 
di(b 23 ) 

d 2 (bi 3 ) 
\ d 3 (b 12 ) ) 



2 * 
-UJ u 



\ A* / 



/ \ 

d 2 H 
d 3 n 
/' 
/' 

fJL 

■dkUi + djUk + d^Uj with k, i,j G {1,2, 3} 
are distinct numbers. Here u* = [u\ + u 2 + u 3 ). 

In the following we will show that how we can get four linearly independent 
vectors of the form v on some subdomain of X. The key thing is to observe the 
behavior of v. We plug the CGO solutions u p = n~ l l 2 w p + fi^Vfp — / p V/x _1 
into v p . Then we get 



(44) 



D ip-x 



( ( -fi- 1 / 2 Pl (p ■ r) - 2/i-Vi(p ■ Va) + 0(|p|) \ 
-/i~ 1/2 p 2 (p • r) - 2^p 2 (p ■ Vs) + 0{\p\) 
-li- l ' 2 P3 {p ■ r) - 2p^p 3 {p ■ Vs) + 0(\p\) 



ip * (p-r) + 2ip 1 p-Vs 
ip^ 1 ^ 2 (p ■ r) + 2ip~ l p ■ Vs 



ip 



-1/2 



(p • r) + 2ip 1 p ■ Vs 



\ 



0(1) 
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Note that r = (n,r 2 ,r 3 ) T and p = (pi,p 2 ,p 3 ) T 

Now we fix any point Xq £ dQ. Let p = r(l,z,0) and p = r(l,0, z) with 

r > 0. Let = p/r, 9 = p/r. Since, in Lemma 2.1, the matrix solutions 

Cq(x,9) and Co(x,9) are invertible, we can choose two constant vectors 

and |(°) such that Co(x, = (r(°),s (0) ) T and 4M)^ (0) = (f {0) ,S {0) ) T 

with sW(x ) = 1 = 5 (0) Oo) and s<°\ |(°) ^ and p • r®, p • f {0) ^ in a 
neighborhood of Xo, say C/q. Then we have the CGO solutions of the elasticity 
system, that is, 

^) = p^(o) +/i -i v /(o)_/(o) v ^i 

with 

Let p = r(i, —1, 0) and 9 = p/r. Let C*i(x, 6 1 ) and C 2 (x, 0) satisfy that 

i9 ■ —C x {x, 9) = 9- V 1 (x)C 1 (x, 9), id ■ — C 2 (x, 9) = 9- ^(x)C 2 (x, 6), 
ox ox 

respectively. Since p = ip, we can choose C 2 (x,9) = Ci(x,9). Moreover, 
r ( 2 ) — r (!) anc i s (2) — s (i)_ With suitable constant vector g, we can get that 
is zero at point xq and r^(xo) = (1, — i, 0) = r^ 2 ^(xo). By continuity of 
r ( fc ) ; we have p ■ ^ in a neighborhood f/i of xo- Then the CGO solution 
is 

u<*> = p" 1 ^ + p^V/f - /f V/i" 1 

with 

For p, with a suitable constant vector g, we can get that is zero at point 
xq and r^(xo) = (1,0, —z). By continuity of r^ 3 \ we have p • ^ in a 
neighborhood £7 2 of xq. Then the CGO solution is 

uf = p^wf + p-V/f - /f Vp" 1 
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with 



w 



(3) 



A3) 

Jp 



D ip-x 



.(3) 
f (3) 



Let U = nl =0 Uk- So far we have five CGO solutions, that is, u p °\ uf\ v!p\u^p\ 



and uf\ 

Using this equality 

(45) 

and (44), we obtain that 



u; 



- %p - x \p\- 2 < 



V 



o 
o 
o 
o 



0(i 



where 



Moreover, at point Xq, 



ul(x ) 



1 y/JI 



2\ + 2fi 



( 2 \ 

2i 





V o / 



0{r- 1 ) 



Similarly, we use (44) and (45), then 



u 2 p := fT 



V^rfirW - r {1) 





0(1) \ 
0(1) 



1 A+2/ t VI 

A V» (JX) 

A / (i) 



2 > 
2 i 

0-h 



and 



Up(^o) 



1 v 7 /^ 
V/2A + 2// 



/ -2r + 0(l) \ 
-2ri + 0{l) 

2/ 
2z 

V 2i ) 



0( 
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Also, we get that 



u 



and 



1/1,2 



1 

71 



u?(*o) 



v 7 ^ 
A+2/.t 

A+2^t 



,(2) 



-zr 



(2) 



(2)> 
2 > 
,(2)> 



1 

\/2 A + 2/x 



' A+2/.t 
' A+2/.t 
' A+2/.t 





^ r ( 2 ) _ r ( 2 )^ 



(2) „(2) 



( 2t + 0(1) \ 
2W + 0(1) 


-2 
-2 

V -2 / 



0(1) \ 
0(1) 



0(- 



.-In 



Let ui = + Since = and s^' = we have 



and 



u, 



1 

71 



/ 



0(1) 
0(1) 


*)rf } + 

2)rf } + 



1 y 7 ^ 

A + 2/^ 



\ 

-l + i)r 2 W ) 
■l + 0r 2 W ) 



0(1) 


2i - 2 
2i - 2 
\2i-2 / 



Oir- 1 ] 



0(1 



The last one is 



u 



5 ._ 



-V- 



f r (3) , „ V (3)n 
A+27ll r l + zr 3 J 



A+2^"i r 3 ; 






Oir- 1 ) 
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with 



1 y/jl 



2\ + 2fi 



( 2 \ 


2i 





V o / 



0{r~ l ). 



Then 5Ju* and Qu p , Jiu^ and $t.u 5 p are linearly independent in the neighborhood 

U of xo as r is sufficiently large. Therefore, there exist functions 9j, Oj and 
j = 1, 2, 3, 4, such that 

3fcr^°/|p| 2 + 9M + 2 9fu* + 6M, + 4 9&u* = 0, 



%e- ip - x v° p /\ P \ 2 + mu p + + mu* + m up = o 



and 



^e-^ x d° p /\p\ 2 + ^SfhiJ + $2%^ + ^u 4 p + = 0, 

respectively 

Summing over j, we have the following equations 
fa ■ V/i + 7iA* 

where 



-U 2 Ui 



for i = l,2, 3, 




-lp-X 



\p\ 



bl 3 
b\ 2 



\ p-ip-x 














-lp-X 



6? 2 



and 7j, Wj are functions which depend on the internal data. 
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Lemma 4.1. /3j is not zero in some neighborhood of Xq, for j = 1,2,3. 

Proof. From above construction, we have knew that 0, ~ r°, 9 ~ r° and 
i?i ~ r°. Applying CGO solutions into the following vector, we get 



P = (P1,P2,P3)- 




With suitable CGO solution u^ ' with r^ ' small enough so that #3(0:0) can be 
as small as we want. Therefore, (3i(xo) is not zero. By continuity of 0i, there 
is a neighborhood V of xq such that /?i never vanishes in V. Similarly, we can 
prove that (3j, j = 2, 3, are also not zero in V. □ 

Let Qo = UdVnQ. Based on the lemma above, we may suppose that fij are 
linearly independent in Q Q as r sufficiently large. Then we have the following 
equation 

(47) fa ■ V/i + 7^/1 = -u 2 Uj, j = 1, 2, 3, 

with {/3j(x)}j = i >2 ,3 a basis in M 3 for every point x G f2o- Since j = 1,2,3, 
are linearly independent in Q , there exists an invertible matrix A = (a^) such 
that (3j = ajkek with inverse of class C°(Q). Thus, we have constructed two 
vector- valued functions T(x), G (C°(f2)) 3 such that (47) can be rewritten 

as 

(48) Vyu + r(x)/i = $(x) in Q - 

Let gp = Up\gQ and gf^ = u P °^\dn be the given boundary data for j = 
0,1,2,3. Let G H 4 (dtt) be the boundary data close to tjf\ g { p j \ 
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respectively, that is, 
ll<7 (i) 



(flu 
9 P \\H*(d<n) 



< e, 



\9 



(o) _ ~(o) 



9 P \\H 4 (dn) 



< e. 



Then we can find solutions of the elasticity system with boundary 

data g( \gti\ respectively. By regularity theorem, it follows that 



(49) \\u^-u^\\ HHQ) <Ce, \\u 

Using Sobolev embedding theorem, we have 



^-uf\\ HHQ) <Ce. 



U)\ 



< Ce, 



(50) u^\\cp(n) 

From (50), for j = 1, 4, 5, we have 



\u 



{3) -~uf\\c^) <Ce. 



$M p \\c°(n) < Ce 



and 



|9fUi - SUp|| c o (n) < Ce. 



Therefore, 9Rui, 9ui, 3?U4, and Jhis are also linearly independent when e is 
sufficiently small. 

Then we have the following uniqueness and stability theorem. 

Theorem 4.2. Let Q be an open bounded subset ofM. 3 with smooth boundary. 
Let (A, fi) and (X,jl) G V. Suppose that [i\qq = jj\gn- Then there is an open 
set of the given boundary data such that the uniqueness and global stability of 
H hold. 

4.2. Global reconstruction of A in 3D. The elasticity system can also be 
written in this form 



(51) 



V ■ u 

V ■ u 

b 23 
bn 

\ 6l2 



difi 
d 2 \i 
d 3 n 



+ 



f d^V-u) \ 
d 2 {V -u) 
d 3 {V-u) 
di(b 23 ) 
d 2 (b 13 ) 
\ d 3 (b 12 ) ) 



( A + \ 

A + fj, 

f 



2 * 
-U U . 



As in the reconstruction of we will construct CGO solutions such that the 
first term of (51) can be eliminated. 



Plug the CGO solution u p = ji 



-1/2 



A* 1 Wp - /pV/i 1 into 



(V-u \ 

V • u 

V • u 

b 23 
613 
hi 
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Then we get 



/ 



D ip-x 



( 



ip l l 2 p ■ r 



ip' 



'1/2 
-1/2 



p ■ r 



2ip~ 1 p ■ Vs \ 
+ 2ip~ 1 p ■ Vs 



p ■ r + lip p • Vs 



0(1) 



-2p~ 1 p l (p l + p 2 + p 3 )s + 0(\p\) 
-2p~ 1 p 2 (p 1 +p 2 + p 3 )s + 0{\p\) 
V V -2/i~V 3 (pi + P2 + Ps)s + 0(\p\) J J 

For the same fixed point Xq E <9f2. We choose a constant vector such 
that C (x, % (0) = s (0) ) with s (0) ^ and p-r^ ^ in a neighborhood of 
x , say U , and p-r^°\x ) = 1. Then we get the CGO solution of the elasticity 
system 



(0) 



- 1 / V°) + p-'vf 



(o) 
p 



with 



(o) 

f (0) 

Jp 



D ip-x 



r (0) 
s (0) 



Of 



We choose another constant vector g^ such that C*i(x, 



(i) 



r 



(1) 



with s^ 1 ) 7^ in a neighborhood of Xq, say Ui, 
get the CGO solution of the elasticity system 



u 



(i) 



-1/2 W (D 



with 



f (i) 
Jp 



D ip-x 



and p 



x ) 



0. Then we 



r (i) 
s (i) 



Of 



For /c = 2,3, we choose a constant vector g^ k ' such that Ck(x,9)g^ = 
(r™, s^) with p-r^ ^ in a neighborhood of xo, say L^- Here we can choose 
r (2) = r ( 3 ) j s (2) = s (3) by taking g (2) = g (3) and C 2 (x,9) = C 3 (x,9). Then we 

get the CGO solution of the elasticity system 



u 



(k) 



with 



+ A* -1 V/« - f { p k) Vp 



(2) 
Wp 

f (2) 

(3) 
Wp 

f (3) 
JP 



^ip-x 



D ip-x 



r (2) 
s (2) 

r (3) 
s (3) 



0( 



0(r 



For p, we choose another constant vector g^ such that C^(x, 9)g^ 
with ^ in a neighborhood of x , say £/ 3 , and p 



r 



(4) „(4)> 



-(4), 



0. Then we 



get the CGO solution of the elasticity system 



u 



(4) 



p- 1/2 w^+p- l Vf^ 



with 
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(4) 
Wp 

(4) 
P 



D ip-x 



-(4) 



0{ 



Then we have that 



u. 







1 + i 
% - 1 




Oir- 1 ) 



For fc = 2, 3, we have that 



u 2 :=e- ip - x \p\- l v 



1 

71 



'A+2u 



p • r 



(2) 



^■r (2) 

-2 /U - 1 r(l + z)s( 2 )+0(l) 
-2 / u- 1 r(z-l)s( 2 ) + 0(l) 
0(1) 



and 



u 3 -.^e-^p^vl 



p ' 



/ 



1 

71 



^p- (3) 

+-p • r (3) 



V 



Vp_~ 

-2p-V(-l +0(1) 
-2/i-V(l -i)s^ +0(1) 
0(1) 



0(i 



Let \x A = u 2 + u 3 then the higher order is eliminated. Thus we have 



u. 



1 

71 



'A+2/u 



0(1) 
0(1) 
0(1) 



+ 0(0. 
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Also, we get that 



4 ■■= e-*np\- 2 vt 



-2// _1 s 



( \ 





l + i 


\i-U 



From above, we have constructed four linear independent vectors Jiu* and 
Qiip, $tUp and JJujj in U = C\\ =0 Uk as r is sufficiently large. Therefore, there 

C>, . 1 I 



are four functions 9 k , k = 1, 2, 3, 4, such that ^e- ip - x v°/\p\ 2 + O^u 1 + 9 2 %u 



8M 



9^u 5 p 



0. Summing over j, we get the following equation 



(52) 



/ X + fi \ 
A + pi 
A + n 

P 
P 
P 



where 



-tp-x 



\P\ 2 

+ m 



-u 



(0)* 



+ 9^ 



-tp-x 



-ip-x 



\p\ 



-U 



(2)* 



IpI 5 



-u 



(1)* 



-tp-x 



9 2 %- 



-tp-x 



\p\' 



-U 



(1> 



\P\ 



-u 



(3> 



+ 04^" 



-ip-x 



\P\' 



-U 



(4)* 



and 



ft- 



-tp-x 



/ ft(v-4 0) ) 

ft(v-4 0) )) 

<9 3 (v-4 0) ) 

ft (68s) 
ft(&? 3 ) 
<9 3 (&? 2 ) 



-tp-x 



/ ft(v-4 1} ) 
ft(v • 

ft (&) 
ft(^ 3 ) 



-tp-x 



\p\ 



( ft(v 
ft(v 

ft(V 



u 



(2), 
P , 
(2), 
P > 
(2), 



ft(6? 3 ) 
ft(&?2 



-tp-x 



+ 



IpI 



) / 



a 



( ft(v 

ft(V 
<9 3 (V 

ft ("2 
ft(&? 3 ) 

V <9 3 (6? 2 ) 



(3)> 
P ; 
(3)> 



-ip-x 



U 



(4), 



9 2 (v-4 4) : 

5 3 (v-4 4) ; 



V 



d 3 (b* 



12 



) / 



ap 



u 2 u, 



We obtain that 
(53) 
where 

k = (i, i, i, o, o, of • (t>° + ^KdJ + 2 3oJ 

a = -(1, 1, 1, 1, 1, if ■ (o° + ^KoJ + 2 3o£ 
Lemma 4.3. k does not vanish in some neighborhood of Xq. 
Proof. Similar argument as Lemma 4.1. Since 



e 3 ^ 

6M 



u, 



«9 3 (V-n p ) 



D ip-x 



w 



and p ■ r^(xo) = 



\+2p 

p • r^ 4 ^(xo), we have 



\+2p 

^p 2 (p-r) 



0(H) \ 
0{\p\) 

o(\ P \) ) 



0(1) 



-ip-x 



-lp-X 



( d 1 (\/-u p 1) ) \ 
d 2 {v- u 

V&(v-t#o y 
/^i(v-4 4) ) \ 



(x ) + # 2 3- 



-tp-x 



/ 9i(V- 
<9 2 (V ■ 

^ 9 3 (v-4 i; ) y 



p 

P 



Oo) 



<9 2 (V-u (4) 



it 



P y 
(4)> 



/ 



\<9 3 (V 

Hence, we obtain that k(xq) is equal to 



(1,1, If 



-lp-X 



( 9i(v-4 0) ) \ 



<9 2 (V • Mp 



(0)> 



v^(v-^ 0) ) y 



+# 3 



-ip-x 



\P\ 

Vp 



( 0!(V- 

<9 2 (V- 

^ 3 (v ■ 4 2) 



(2)^ 



-ip-x 



(x ) 



+ 



)y 



<9 2 (V-« 

V ^(v-«n yy 



(3)> 
P > 
(3)> 



(x )) + O(r- 



/ 



2(A + 2p) 



-lp-X 



\p\ 



-tp-X 



i 0i(V- 



+ 



(3)> 



<9 2 (V ■ u (3) 



\ <9 3 (V 



P ; 
(3)> 



yy 
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With suitable chosen CGO solution with small, then it implies that 
#3(^0) can be as small as we want. Thus, k(xq) 7^ 0. By continuity of k, k 
does not vanish in some neighborhood V of xq. □ 

Let Qo = U fl V n Q. Then k does not vanish in Qq. Then it follows that 
(54) A = ^ - to 2 - in fi . 

Then we have the following theorem. 

Theorem 4.4. Let Q be an open bounded domain ofM? with smooth boundary. 
Suppose that the Lame parameters (A, fj) and (A,/i) G V and /i|an = Alan- 
T/ien i/iere zs an open se£ of the given boundary data such that the uniqueness 
and the global stability of A hold. 
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